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Abstract 



m 

■ We analyse qualitative properties of the solutions to a mean-field equation for particles interacting 

00 , through a pairwisc potential while diffusing by Brownian motion. Interaction and diffusion compete 

OO ' with each other depending on the character of the potential. We provide sufficient conditions on the 

relation between the interaction potential and the initial data for diffusion to be the dominant term. 
We give decay rates of Sobolev norms showing that asymptotically for large times the behavior is then 
given by the heat equation. Moreover, we show an optimal rate of convergence in the L^-norm towards 
the fundamental solution of the heat equation. 



^ I Mathematics subject classification numbers: 35B35, 35Q30, 76D05 

1 Introduction 

In this paper we consider the diffusive aggregation equations, 

dtp = V ■ {p{VW * p)) + Ap (1.1a) 
p{0,x) = po, (1.1b) 



where p = p{t,x) is a real function depending on time t > and space x G M^, W : — > M is an 
interaction potential verifying W{x) = W{—x) (without loss of generality, see [15]). These equations 
have received a lot of attention in the recent years because of their ubiquity in different models and areas 
of applied and pure mathematics. Collective behavior of animals (swarming), chemotaxis models, and 
granular media models are some examples, see [26, 27, 11, 6, 15] and the references therein. On the other 
hand, these equations have been studied in connection to entropy-entropy dissipation techniques, optimal 
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transport, and gradient flows with respect to probability measure distances, see [15, 2] and the references 
therein. 

Without diffusion, the continuity equation (1.1a) with a singular interaction potential W can lead to 
very involved dynamics where blow-up can occur, and where Dirac Delta singularities and smooth parts 
of the solution can coexist. More precisely, assume that we have an interaction potential which is radial, 
smooth away from the origin and whose gradient may be singular at the origin, with a local behavior not 
worse than that of |x| (i.e., having at worst a Lipschitz singularity at the origin); then, blow-up in finite 
time of n solutions was reported in [7]. In fact, the almost sharp condition which determines the 
behavior of global existence or blow-up in n is the so-called Osgood condition. It is given in terms 
of the size of the gradient of the potential W, which is said to satisfy the Osgood condition if 



with W{x) = k(\x\). Specifically if (1.2) is satisfied, with some mild additional monotonicity conditions 
on k", and also VW £ W'^'" with q < N and po e D LP with p > N/{N - 1), N > 2, then there 
are global weak solutions, see [8]. Aggregation (blow-up) only happens when t = oo, see [7, 8] for the 
L^-L°° and L^-U' results, respectively. In fact, weak measure solutions were proved to exist after the 
L°°-blow-up time in a unique way for certain attractive non-Osgood potentials, see [13] for a definition 
of weak measure solution and further details. In [13], the authors also illustrate the existence of weak 
measure solutions with very complicated patterns and Dirac Delta formations. 

In this work we address the following issue: under which conditions on the interaction potential can 
linear diffusion prevail, leading to a diffusive-dominated behavior for large times? More precisely, we give 
sufficient conditions on the interaction potential and the initial data such that the competition between 
the possible aggregation due to an attractive interaction potential and the linear diffusion in (1.1a) is 
won by the latter. 

Sharp conditions for separating global existence of solutions from blow-up of solutions to (1.1) have 
been given in several papers related to homogeneous interaction potentials or to the classical Keller-Segel 
model [11, 9, 5]. Also, conditions for global existence or blow-up have been given in [22] based on the 
regularity of the gradient VW of the potential. Blow-up conditions have been studied in detail as well 
for fractional diffusions [24, 23]. However, we are not aware of many results dealing with the asymptotic 
behavior once the diffusion dominates over the aggregation except for [10, 4] where the authors show that 
the solutions of the Keller-Segel model behave like the solution of the heat equation for small mass and 
the recent papers [21, 22] in which they deal precisely with this issue for the equation (1.1); we comment 
further on these works below. 

Here we show that under suitable smallness conditions involving both the interaction potential and 
the initial data, see Theorem 3.1, the behavior of the solution is determined by the heat equation for 
large times (see Theorems 3.1, 3.5 and 4.1). In other words, we get a result of asymptotic simplification 
for all dimensions under some size conditions where the nonlinearity disappears and the decay rates 
and behavior are like the diffusive equation at least up to first order. We remark that the asymptotic 
simplification result in without rate and the decay rates in were obtained in the one-dimensional 
case in [21] under some smallness condition similar to one of the possibilities in Theorem 3.1 below by 
using scaling arguments. Also, global existence results were reported in the multidimensional case in [22] 
but no uniform-in-time bounds nor decay rates under smallness size conditions were proved. 




(1.2) 
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Asymptotic simplification results have been reported for problems in fluid mechanics [1, 28, 29, 20], 
in convection-diffusion equations [14], and in some nonlinear diffusion models of Keller-Segel type for 
small mass [4, 25]. Here, we use the technique of Fourier splitting [30], a technique quite successful for 
the > 3 dimensional Navier-Stokes equations, together with direct estimates over the bilinear integral 
term associated to (1.1a) via Duhamel's formula to get the optimal decay rates in Sobolev and spaces 
in section 3. Section 2 is devoted to setting the basic well-posedness theory of global- in-time solutions 
with uniform-in-time L°° estimates. Finally, Section 4 is devoted to combining these time decay estimates 
with entropy-entropy dissipation arguments [31, 33, 3, 16, 17, 14, 15, 4] to obtain decay rates in entropy 
in self-similar variables and in L} towards the self-similar heat kernel for large times. 



2 Well-posedness and global bounds 
2.1 Notation and preliminaries 

We usually omit the variables of the unknown p in eq. (1.1), which are understood to be {t,x). Also, we 
usually write pt{x) = p{t,x), which is useful when referring to the function x i— )• p{t,x) for a given time 
t (we emphasize that pt is not to be confused with dtp; subindex notation for partial derivatives is never 
used in this paper). When not specified, integrals are over all of M^, and in the variable x. 

We use the standard multi-index notation for derivatives throughout: for a function / : M.^ — )• C and 
a multi-index 7 = (71,72, ...,7Ar), with integers jj > 0, we denote I7I = X^^i 7j and define 

B"' f = f)^^ rTi'^ Ff<^ f 

We let S = 5(M^) be the usual Schwartz space of rapidly decreasing functions. The Fourier and inverse 
Fourier transform of / G 5 are defined by 

/(^) = (2^)-^/2 / e-'^'-^vix) dx and f{x) = (2^)-^/^ / e^^-«/(0 d^, 



respectively, and extended as usual to S'. If /c is a nonnegative integer, W'^'^(M^) = W'^'^ will signify, 
as is standard, the Sobolev space consisting of functions in L^(M^) whose generalized derivatives up to 
order k belong to LP(R^) = L^, I < p < 00, with norm || • [[^p defined by 

wfwip ■■= E w^yWp (2.1) 

7I <fc 

for any / E VV*^'^, where the sum is over all multi-indices 7 with I7I < k. When p = 2, W^'^(M''^) = 
H''{R^) = H^, where the space H'^ is defined for all /c € M as the space of all f € S' such that 
(1 + |C|^)'^/^/(0 £ L"^- The norm in H'^ (that of W'^'^ defined in (2.1)) is sometimes denoted by || • W^k 
instead of || • ||fc^2- We also use the notation 

The space of bounded continuous functions from an interval / to a normed space X is denoted by 
BC{I,X). To simplify the notation, the domain of the above spaces will usually be omitted. 
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We will use the well-known fact that for any m > there is a constant C depending only on m and 
the dimension N such that 

\\f9\\m,p < C{\\f\\mj9\\oc + ||/||oo||5lU,p), (2.3) 

for functions f,gG W"*'^ H L°°. We also define D^-p by D'^p = \^\'^p. Hence, by Plancherel's identity 
||-D"^/>||2 = |||CI™'p||2. By C we denote arbitrary constants that can change from line to line. 

We recall the Gagliardo-Nirenberg-Sobolev (GNS) inequalities, which we will be using repeatedly in 
the sequel, see [19, Theorem 9.3] for a proof: given 1 < q, s < oo and integers < j < m, there exists a 
number C > depending on q, s, j, m and the dimension N such that 

\\V^v\\p<C\\V'"'v\\l\\v\\l~^. (2.4) 

where we use the notation (2.2) with 

1 ?' „/l m\ , 1 

r^+''(^-«)+<'-'''? 

where ^ < < 1, with the following exception: if m — j — N/q is a nonnegative integer, then the 
GNS inequality (2.4) is only valid for ^ < < 1. Given that the seminorms ||Z)™u||p and HD"*?; ||p are 
equivalent for 1 < p < +oo, we will often interchange them. 

We will make use of the standard heat semigroup e*^, which is defined as the convolution in the x 
variable with the heat kernel ^ 

G(^,x):=^^e-^, (2.5) 

whose derivatives satisfy for all integers m > 1 and any multi- index 7 with I7I = m 

_/V"4-m \x\'^ 

\d"'G{t,x)\ <Ct — —e~^ (2.6) 



for some constant C > depending on the dimension A^. From (2.6), we get the following standard 
estimates: 

\\VG{t, ■)\\r < Ct-i(i+^(i-^)) (2.7) 
with 6* = i and C = C{N) and 1 < r < oo. 



2.2 Local-in-time existence 

In order to show short-time existence of solutions to the equation (1.1) we use a common fixed-point 
iteration as also done in [22]. Some elements in the proofs of subsections 2.2 and 2.3 are related to results 
reported in [22] but we prefer to include them for the sake of the reader. One can formally rewrite (1.1) 
by using Duhamel's formula and integrating by parts: 

Pt = e*^po - r Ve(*-^)^ {psiVW * ps)) ds, (2.8) 
Jo 

where Ve*^ denotes the convolution in x with the gradient of the heat kernel (2.5). As it is common, we 
define a mild solution of (1.1) as one that has some reasonable regularity for (2.8) to make sense, and 
satisfies (2.8). For this definition, and for the short-time existence result 2.4, we do not assume that p is 
nonnegative, as it is not needed for the argument. In the sequel we work in any dimension > 1. 
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Definition 2.1. Take T G (0,+oo], p G [l,+oo], and po G L^. Assume that VW G (L^)^. A mild LP 
solution to equation (1.1) on [0, T) wit/i initial condition po a function p G LjQ^([0, T), L*') stic/i i/iai 
(2.8) holds for all t G (0,r). 

Remark 2.2. iVote i/iai VW * p G L/^^([0, T), (L°°)^), so if p e LP the product p{VW * p) is in 
G L;^Q^([0, T), (L^)^) and the integral in (2.8) makes sense. 

Remark 2.3. is easy to see that whenever a mild solution p has enough regularity to he a classical 
solution, it is then in fact a classical solution. This is: if a mild solution p has continuous first-order 
time derivatives and continuous second- order space derivatives, then it is a classical solution. 

We will show the following short-time existence theorem: 

Theorem 2.4 (Short-time existence). Take p G [1,-|-cxd], m > 0, and po G W^^'P. Assume that VW G 
{L'^)^ , with 1 + 1 = 1- Then there exists a maximal time T* G (0, -|-oo] and a unique mild solution 
p G C([0,r),W™'P) of problem (1.1). If T* < +oo then 

\\pt\\m,p ^+oo as t —?■ T * . (2.9) 

Let p' G [0, +co], m' > 0. // additionally po G ^^''P' and VW G (L"')^ with + ^ = I, then the 
solution given above belongs to BC{[0,T),W"'^p HW^'^p'). 

This result will follow from a standard fixed point theorem for bilinear forms which for completeness 
we state here, see [12]: 

Lemma 2.5. Let X be an abstract Banach space with norm \\-\\x oind B : X x X ^ X a bilinear operator 
such that for any xi,X2 G X, 

||-B(xi,X2)||x < ??||xi||x||x2||x (2.10) 

then for any y £ X such that 

Ar]\\y\\x < 1 (2.11) 

the equation x = y + B{x,x) has a solution x £ X. In particular the solution satisfies \\x\\x < 2||?/||x 
and is the only one such that \\x\\x < 

Proof of Theorem 2.4- Take T > 0. Following a standard strategy, we first show that for T small enough 
there exists a mild solution on [0,T). In order to find a function p satisfying (2.8) we choose the bilinear 
form defined as ^ 

B{p, 1^) = - [ Ve^^-''^^ ■ {VW * il)s)ps ds for t G [0, T). (2.12) 

In order to apply Lemma 2.5 estimate this bilinear form in the space X := BC{[Q, T), W™'*'), with norm 
given by 

\\p\\m,p ■= sup \\pt\\m,p- 

t€[0,T) 
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In this proof we denote by Cm a number that depends only on m, which may change from hne to hne. If 
7 is a multi-index with I7I < m, then, for p,ip £ BC{[0,T),yV^'P), and any t G [0, T), using (2.7), we get 







\\d"'B{p,^l;)\\p< I Ve(*-^)^-9^((VVF*V'.)Ps) < C j -j== 5t((VVF * Vs)p.) 



p Jo 



* 1 



ds 

p 



/■* 1 „ 

<Crn —^^=\\VW\\q\\'tps\\m,p\\Ps\\m,pds, (2.13) 

Jo Vt- s 

which shows, taking the supremum on [0, T) and summing over all multi-indices with I7I < m, that 

mp,mm,p < CmVT\\vw\u\Mm,p Us\im,p- 

In the third inequality in (2.13) we have used that, for / G W™'^, g G 

W{f9)\\p<Cm\\f\\m,p\\9\U,oo, (2.14) 

for any multi-index with I7I < m. This can be easily seen by writing out and estimating the derivatives 
of the product fg. This gives the estimate (2.10) with r] = CmVT \\VW\\g. Taking y G BC{[0, T), W™'^) 
defined by t 1— )• e^^po in Lemma 2.5, we can choose T small enough so that (2.11) is satisfied. Hence 
Lemma 2.5, yields the existence of a function p G BC{[0,T),W^'^) satisfying (2.8), i.e., a mild solution. 
This solution is a priori unique only in the set of solutions satisfying < 1/(2?/), but a standard 

argument using the continuity of p shows that it is in fact the unique mild solution in BC{[0,T),W^'^). 

The existence of a maximal time T* and the blow-up of the solution at T* if T* < +00 follows 
now from a standard argument. The last part of Theorem 2.4 is obtained by an analogous reasoning, 
considering now the space X := BC{[0, T),W"''P n W^'-^'). □ 

2.3 A priori time-dependent bounds and global existence 

In this section we obtain bounds for our solutions. Due to equation (2.9) in Theorem 2.4 these bounds 
imply global existence of the solutions. 

Proposition 2.6. 

i) Letpo eL^riLP with po > 0, andVW G (L^nL^)^ with l/p+l/q = 1. Let p C([0, T); n L^) 
he a mild solution to (1.1) on [0, T) as obtained in Theorem 2.4 with initial data pQ. Then there is 
a constant C > ( depending only on pQ and the dimension N ) such that for all t G [0, T) 

\\pt\\i = \\po\\i=- M, (2.15) 
WptWp < ll/oollpexp (c\\VW\\oo\\po\\iVi) := Cp{t) ||po||p- (2.16) 



a) Assume that po G W"'^ D D L°° for some m > 1, and VW G {L^ D L°°)^ . Then there is a 
time- dependent function Cm,p{t), hounded on finite time intervals, depending only on po, ||VVF||i, 
||VM^||oo and the dimension N, such that for all t G [0,T), 

llptll/f™ < Cm,p{t) IIpoIIh™- (2-17) 
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To prove this proposition we will use the following modified Gronwall Lemma (see [29]): 
Lemma 2.7. Let <T < oo, 6 £ (0, 1), and let f : [0, T] — )• [0, oo) be continuous and satisfy 

f{t)<A + B [\t-sr'f{s) ds 
Jo 

for all te[0,T). Then 

f{t) < A<!>{BT{1 - 6)t^~^) 
for t E [0, T) , where <^ : C — )• C is defined by 



Proof of Proposition 2.6. Inequality (2.15) follows by direct integration of equation (1.1) (or of the mild 
formulation (2.8)) since it has divergence form. To obtain (2.16), take the norm in (2.8) and use (2.7) 
to obtain 



(2.18) 



/■* 1 

Ilptllp < IIpoIIp + C" / —==\\ps\\p\\VW * ps\\oods 
Jo Vt- s 

/■* 1 

<||po||p + C' / WpsWp ||VT4^||oo||p.s||i ds 

Jo vt — s 

„ „ r* 1 „ „ 

< IIpoIIp + C'llVWIIooll/Oolli / -i^=\\Ps\\pds. 

Jo Vt- s 

The modified Gronwall inequality (Lemma 2.7) yields (2.16). To obtain (2.17) take any multi-index 7 
with I7I < m, apply d'^ to (1.1) and multiply by d'^p to obtain 

(To make this reasoning rigorous, as we are using eq. (1.1) instead of the weak formulation, we have to 
carry it out on approximating solutions with smooth initial data and then pass to the limit. This process 
is straightforward and as such we omit the details.) By Cauchy-Schwarz's inequality the first term can 
be bounded by 

J V{d'rp).{d'r{p{VW*p))) < (I |V5VP)'^'(| \d^{p{VW*p))\'y\ (2.20) 
To estimate the second parentheses, use equation (2.3) and then Young's convolution inequality to obtain 

/ \d^p{VW * p))n < \\p{VW * p)\\m,2 < C{\\p\\m,2\\'^W * p||oo + ||/o||oo||VVF * p\\m,2) 

< C\\VW\\^ ||p|U,,2(||/o||i + IIpIIoo) =: c{t)\\ 

p\\in,2i 

with C{t) db givGn function that involves ||VW^||oo; IIpoIIi and Coo{t) from (2.16). Using this in (2.20) and 
applying Young's inequality we get 

I V(a» • {d^{p{VW * p))) <\j |VaV|2 + C{tf\\p\\l^^. (2.21) 
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Combining this with in (2.19) yields 
Adding all multi- indices 7 with I7I < m gives 

J^\\p\\l,2<C{t)^pC^^ 

for some other time-dependent function C{t). Integrating this inequality over time proves the last part 
of the proposition if all derivatives above are well defined. □ 

Now, we combine the short time existence in Theorem 2.4 with the a priori results in Proposition 2.6 to 
yield the global existence. 

Theorem 2.8. Under the conditions i) of Proposition 2.6, there exists a unique global mild solution p of 
(1.1) with p £ C([0, 00); LP). Under the conditions ii) of Proposition 2.6, there exists a unique global mild 
solution p of (1.1) with p £ C([0, 00); ff"*). 

2.4 Uniform bound for ||p||oo 

Consider the solution p G C([0, +00), n L°°) obtained from Theorem 2.8 with p = 00. In this section 
we prove that the L°° norm is actually uniformly bounded for all times: 

Theorem 2.9. Let the interaction potential W be such that VW G (L^ n L°°)^ . Let po £ L^ D L°° 
nonnegative, and suppose p is the solution constructed in Theorem 2.8 with data pQ. Then there exists a 
constant Coo depending only on N , W , and ||po||i such that 

llPtlloo < Coo for all t > 0. (2.22) 

Proof. Choose any time to > 0. From Proposition 2.6 there is a constant Cq for which 

||aI|oo<C'o, tG[0,to]. (2.23) 

We will prove that, for some 6 > 0, 

llPtlloo < Coo , t£ [to -6,00). 

This will ensure the global-in-time bound for the L°° norm. 

Pick 5 satisfying to > ^ > (further conditions on 6 will be fixed below). As p is a mild solution, we 
may use Duhamel's formula (2.8) between t — 5 and t to obtain, for any t > 6, 

pt = e'^pt-s - f Ve(*^^)^(p,(VVF * ps)) ds. (2.24) 
Jt-& 

For the first term in (2.24) we have 

e^^P*-. < / e-^p(t - 5, y) dy (2.25) 

1 1 

- (47r(^)^/2"^*-'^"^ " (47r<5)^/2"^°"'- 
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To bound the second term in (2.24) let g := {VW * Ps)ps-, 

[ Ve(*-^)^-5ds < / \\VG{t-s,-)*9\\oods< I \\VG{t - sr)\\r\\9\\q ds (2.26) 
Jt~5 Jt-5 Jt-S 

with i + i = 1. We use the fohowing estimate for H^Hg: 

\\g\\q < \\Ps\\oo\\VW * 

PsWg ~ \\Ps\\oo ll'^^llg IIpoIIi ! (2.27) 
and (2.7) for \\VG{t - s, OUr with 6^ = i, to conclude that 

f Ve^'-^^'^-gds <C||VT^||,||po||i t (i - .)"^(^+^(^-''» ||p.||oo . 
Jt-S Jt~5 

To finish the argument it is necessary that the right-hand side of the last expression is integrable near 0, 

hence we need 

1 1 

_ 1 < __(1 + iV(i _ 0)) ^ 1 _ < _ ^ ^ < (2.28) 

which means, from (2.28), that we need N < q < oo. Putting (2.25), (2.26), (2.27) and (2.7) together in 
(2.24) we obtain, for t > 6, 

\\pt\\oo< ,^ l,^M i + c\\vw\u\\poh f (t-s)-^('+^('-'))||Ps||oods 

{ATTdy^l'^ Jt-5 

<Cr^/2||po||i + CC2 sup \\p{r)U f [t-s)--^^^+^^^-<^)Us 

Te{t-S,t) Jt-S 

<Cr^/2||po||i + CC2 5'^ sup ||p(t)||oo 

Te{t-5,t] 

where C2 '■= ||VVF||g ||/0o||i, C is a generic constant depending only on and 2p = 1 — N{1 — 9) > 0. 
Now, we take the supremum of both sides of the inequality for t > to to obtain 

sup\\pt\\oo<C6-''/^\\po\\i + CC2 6^' sup ||p(r)|U 

t>to t>to-S 

< C <5-^/2 ii^^ii^ + C7 C2 <5^ Co + C C2 5^^ sup ||p(t) lU 

t>tQ 

by using the local estimate in time (2.23) for times less than to- Choosing < < to such that 

2 C C2 = 1 we have 

i>io 

This completes the proof of the theorem. □ 

3 Algebraic decay of the solution 
3.1 decay 

To obtain the decay in of the solutions, we first need to show that solutions satisfy an estimate of the 
form 

J dx < -C J I Vpp dx, (3.1) 

for some C > 0. This will hold for sufficiently small potentials. From this, given that the solution remains 
in L^, decay at the same rate as solutions to the heat equation will follow. 
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Theorem 3.1. Let the potential W be such that VW G (L^ n and take po G n n 

nonnegative. Assume additionally that one of the following smallness conditions holds: 

N+i 

i) It holds that || VH^||oo||po|ir^'' < 1- 

a) W £ and ||H^||iCoo < 1, where Coo such that (2.22) holds (i.e., ||pt||oo < Coo for all t>0). 
Hi) W £ L'^ and \\W\\2C2 < 1, where C2 is such that \\pt\\2 ^ C2 for all t > 0. 
iv) C||po||i||[^W^]+||Af/2 ^ 1/4 where C is given below in the proof and [AW]+ := max{AM^, 0}. 
Then there exists a number /C > depending only on the constants defined in the hypotheses such that 

\\pth<JC{t + ir''/\ (3.2) 

1 1 1 /2 1/2 

Remark 3.2. The constant C2 can be estimated by interpolation as C2 < llpolli Coo ■ However, the 
constant Coo depends both on W and po in a complicated way, although an explicit bound can be extracted 
from the proof of Theorem 2.9. 



d 
di 



Proof. We first establish (3.1) in all cases. Multiply equation (1.1a) by p and integrate to get 

■ j p^dx = j ■ (PV^ * P)'ix - j \^P?dx. (3.3) 

We need to bound the first term on the right-hand side of the last equation. 
Case i): We use Holder's and Young's inequalities to obtain 

j pV ■ {pVW * p)dx = - j Vp- {pVW * p)dx < \\Vp\\2\\p{VW * p)\\2 

< llVplb ||p||2 \\yW * pIIoo < ||Vp||2 ||p||2 llVt^lloo \\p\\l 

JV+4 

<l|vHlil|vw||oo||polir% 

where we interpolated the norm of p between and by means of the GNS inequality (2.4) with 
j = 0, p = q = 2 and m = s = 1 and the conservation of mass (2.15). Combining this estimate with (3.3) 
yields (3.1), due to the condition i) in the hypotheses. 
Case ii): By a similar reasoning we have 

j pV ■ {pVW * p)dx = - j Vp- {pVW * p)dx 

< ||Vp||2||/0(T^*Vp)||2 < ||Vp||2||p||oo||l^||l < \\Vpf2 Coo\\W\\l. 

Combining this estimate with (3.3) yields (3.1), since Coo ||^||i < 1- 
Case ill): Similarly, 



J pV ■ (pVW * p)dx = - J Vp - {pVW * p)da 



< ||Vp||2||p(i^* Vp)||2 < ||Vp||2||p||2||l^* Vpllc 

< ||Vp||2 C2IIVFII2. 
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This gives (3.1) as before, since C2 ||VF||2 < 1 by hypothesis. 

Case iv): Since the Laplacian of the interaction potential Ues in an space, we can write 

j pV • {pVW * p)dx = - j pVp ■ {VW * p)dx = j ^(ATy * p)dx < j y ([AVF]+ * p) =-. I . 

Case iv.a): We first consider the case N = 2. Then, choosing any 1 < p < +00, 

2/ < ||/)||i||[AW^]+ *p|L with - + - = 1. 
' P Q 



By the GNS inequahty (2.4) we have 

\\p\\lp<C\\Vp\\np\\f'~''\ (3.4) 
where a = Now we estimate the convolution by 

\\[AWU*p\U<\\[AWU\\i\\ph- (3.5) 
Interpolating L'^ between and by the GNS inequality (2.4) with b = ^ yields the estimate 

llplk/<C||Vp||^||p||^^ (3.6) 

Note that a + | = 1. Thus combining (3.4), (3.5), and (3.6) gives 

I<2C\\[AW]+h\\Vp\\l\\ph. 

Note that p can be chosen to be any p G (l,oo), and this final inequality does not depend on the choice 
(except for the constant C in front of the inequality). This concludes the proof of the case N = 2. 

Case iv.b): We now consider > 3. In this case the nonlinear term is bounded by 

1 4 2N 

^< 2ii/'iiiii[^^]+*/'ii°o<c^iiHir'iivpiir'ii[Aw^]+*piioo, 

by means of the GNS inequality (2.4). To estimate the term ||[AM^]_|_ * p||oo) we use Holder's inequality 

\\[Awu*p\\oo<mwu\\p\\p\\g, 

with 1 < p < +00 and | + | = 1. We interpolate the L'^ norm between and using again the GNS 
inequality (2.4) to obtain 

llp||.<l|Vp||5||p||}-", (3.7) 

with I = a(^ — j^) + 1 — a, and thus a = p(^^2) • ^ 1 we need p > We want to choose p so 

that 

N a _ N ^ _ 

2 + N 2 ~ JTn ^ p{N + 2) ~ ^' 

This will hold if p = N/2, N > 3. Notice that the GNS inequality (3.7) does not hold for N = 2 since we 
should take a = 1 and q = 00, which is not allowed. From the above inequalities, it follows that 

I<C\\pMVp\\l\\[AW] + U,2. 
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Hence, Case iv) for > 3 follows. 

We now finish the proof of the L^-decay in the three cases above. The GNS inequality (2.4) once more 
gives 



JV 2 

lPll2<l|vp||^^lHlr^ 



JV+2 



Thus from (3.1) it follows that 



d 



1 



— / p^dx < — 
dt ' ^ - 



so ||Vp||2> ||p||2^ 



JV+2 

\Vp\'^dx < -JC\\p\\2'' 



\p\\i 



Integration of this differential inequality yields the expected decay. 



□ 



Remark 3.3. Some typical examples of interaction potentials in applications are variations of the so- 
called Morse potential W{x) = 1 — e"'^'" with a > 1, see [7]. For instance, we can get decay for solutions 
to the aggregation equation with W{x) = 1 — e"'^' provided 



c<^{\\poU[Awuy/,)-' 



here \\[AW]. 



,N/2 
\n/2 



e 2 



" ' ' dx. 



3.2 H"" decay 



In this subsection we consider the decay in H"^ spaces. The aim is to show that for certain potentials the 
decay rate will be the same as the one for the solutions to the heat equation. Specifically, we show that 
if the decay happens at the same rate as for the heat equation (as was shown in Theorem 3.1 under 
some additional conditions) then the decay in H"^ will happen at the same rate as for the heat equation. 
Our potential will satisfy the hypotheses given in Theorems 2.8 and 2.9, which ensures the existence of 
a unique global solution, L°°-uniformly bounded in time. In this subsection, we remind the reader that 
D"^p is defined via the Fourier transform by D"^p = \S,\'^p, as remarked in Section 2.1. 

We first give a technical lemma to be used later in the proof: 



Lemma 3.4. Take m > 1, and functions f,g,h G such that f and D"^g are in Lp- 
D'^{f{g * h)) is also in and 

\\D^{f{g * h))h < 2^-' (||I?'"/||2||<7ll2||/i||2 + \\fh\\D"'9h\\hh) . 
Proof Denote u := g * h. Using that < 2^-1(1^ - ??p + |?/p) for any C,u£C, 



Then 



\D"'ifu)\^dx= / \\^rfu\^dC 



'f*ufd^ 



\crfii^m-i^)d,. 



d^ 



< 2 



2m-2 



l2m-2 



i^r\fi'^)\m-'^)\d>^+ / \ fim-'^rm^-'^)\di^ 



dc 



di 



(ien/i)*N + i/i*(i^rN) 



Consequently, using Young's inequality for convolutions, 

\\D'^{fu)h < 2^-^ (\\{\ir\f\) * \u\\\, + \\\f\* iicnuDh 



< 2' 



■m— 1 



(Iii?'"/ii2iiniii + ii/ii2iiier2i 
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The proof is now finished by noticing that 

ll^^lli = WdHi < Il?ll2||/i||2 = lls-lbll/ilb 

and similarly 

mruWi = mrdhWi < mrduhh = ii^"^<7ii2iiuii2. 

□ 

Theorem 3.5. Let the potential W satisfy VW G (L^ n L°°)^ and let po e D L°° D H"^ , m>l, with 
Pq nonnegative. Consider p the solution to (1.1) given by Theorems 2.8 and 2.9 with data po . Assume 
that the solution satisfies the -decay estimate (3.2). Then there exists a constant C > which depends 
only on W , po, m and N such that for all t > 

WD'^pth < C{t + i)-W4+"V2) , (3_8) 

Proof. In this proof we denote by C any nonnegative number that depends only on the same quantities 
as the constant C in the statement. We first need to show that for t large enough 

^ J \D'^p\^dx <-\j \D'^^^p\^dx. (3.9) 

For this, apply the operator D"^ to the equation (1.1a), multiply by D"^p, and integrate in space. After 
reordering and integration by parts it follows that 

d 



dt 



I \D'^p\^dx = -2 I D'^pD'-^V ■ip(yW*p))dx-2 j ^IdjD'^pl^dx . 



Integrating by parts in the first integral on the right-hand side yields 

d r r ^ r ^ 

— / \D"'p\'^dx = 2 / ^{djD"'p) D'^ipidjW * p))dx - 2 / ^\djD"'p\'^dx. 
•' i=i •' i=i 

By Holder and Young's inequalities, we obtain 

^ j \D"^p\'^dx< j \D"'{p{W*Vp))\'^dx- j Y,\djD"'p\'^dx. 



It follows that 



^J\D'^p\^dx<J\D'^{p{W*Vp))\^dx-J iD^+VPdx, (3.10) 



where we have used that 



N \ / N \ N 



,i=l / \j=l / k=l 



We need to bound the first term on the right-hand side of the inequality (3.10). Using Lemma 3.4, 

I \D^'\p{W * Vp))|2 dx<C {\\D"^p\\l\\W\\l\\Vp\\l + \\p\\l\\W\\l\\D^-^'p\\l) =: / + II. (3.11) 
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To estimate /, we proceed by means of the GNS inequality (2.4) to get 

2m 2 2 2m 

< C\\D"^+'^p\\^+' \\p\\^+' and llVplli < Cp'^+Vlla'"^' llpllr^' 

which yields 

I<C\\D^-^'p\\l\\p\\l\\W\\l 
Since the term II is already of the same type, from (3.11) we then have 

J |z)™(p(VF * vp))p dx < cp'^+Vlli IIpIII IIH^IIi 

Since W ^ and ||p||2 — due to Theorem 3.1, it follows that there exists Tg so that for all t > To, 

j \D"'[p{W * Vp))|2 dx < 1/2||D'"+V||i ■ 

Using this in (3.10), inequality (3.9) follows for t > T = max{To,Ti}. 

Now, after taking the Fourier transform and applying Plancherel's Theorem it follows from (3.9) that 



1 1 \\^rp\'d^<-^ I w^n'pi'd^. 



Now we proceed by Fourier splitting, see [30]. We split the frequency domain into 5 and S'^, where 

2k 



s{t) = u : lei < Gm, Q{t) 



+ 

and /c is a positive number to be chosen later. By Plancherel's identity we obtain 

^ f\D^p\^di<-\ f \i\^\D^p\'^di<--^ f i^pde 



<--^[ \D^p\'dC + -^l \D"^p\'d^. 



S(t) 



Hence, using Theorem 3.1, 

\t + l)''\\D^p\\l\ <A:(t + l)'=-i / W^rpl^d^ 
L J Js(t) 



d_ 

It 



Sit) 

< k (2^)"^ {t + i)fc-i-'» I < cwpoWi {t + pfc-i-'" 

Jsu) 



\r I _ - iir VII i \- ■ -/ ^ 

S{t)"' 

Choosing now k > N/2 + m and integrating on [T,t] gives the desired decay rate (3.8). Together with 
the a priori estimate at time t = T obtained in Proposition 2.6, this concludes the proof. □ 

3.3 LP decay 

We give first the decay for the norm. The decay of the solutions in all other spaces follows by 
interpolation. 
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Lemma 3.6. Let the interaction potential W be such that VW G {L^nL°°)^. Let po e n n 
with m > N /2, and po nonnegative. Consider p the solution to (l.la)-(l.lb) with data po given by The- 
orem 2.8 with the properties in Theorems 2.9 and 3.5 under one of the additional smallness assumptions 
in Theorem 3.1. Then there is some constant C > such that for all t > and p G [2, oo], 

||Allp<C(t + l)-^/2(i-Vp). (3.12) 

Proof. Let us start by the case p = oo. Using the GNS inequality (2.4) with j = 0, a = oo, 6 = 2 and 
s = 1, we obtain 



WPWoo ^ ^ \\^ P\\2 WPWl ■ 

vahd for any m > N/2. Using the decay, we have for ||Z)™73||2, 



IpIIoo < c{i + ty^-+-'^^ = c(i + 1)-^. 



The general case p G [2,oo) follows by interpolating L^ between L and L°°. □ 
3.4 Decay of ||a;p||2 

We need to study the behavior of the norm ||x/o||2, as this norm will appear later in estimates. We begin 
by studying a moment of p: 

Lemma 3.7. Assume the conditions from Lemma 3.6. In addition, suppose that \x\^pQ € L^ and that 
for some C > 0, 

\x-VW{x)\<C , xGM^. (3.13) 
Then, there is a constant C = C(||po||i)-^) such that for allt>d 

j \x\^p{t,x)dx < j \x\^ po{x) dx + Ct and \\xpt\\l < C{1 + tf'^ . (3.14) 

Proof. Multiplying (1.1a) by |xp and integrating, we get 

j \x\'^p{x) dx = 2iV||po||i - 2 J p{x- {VW * p)) dx 

= 2N\\p4i - jjix- y)p{x)p{y)VW{x - y) dxdy < 2N\\p\\i + C\\p4l, 

which establishes (3.14). Now, using (3.12) and (3.14), we get 



\xp\\l 



\x\'^p{xf dx < IIpIIoo / \x\'^p{x) dx < C(l + tf~^ . 



□ 



4 Asymptotic simplification towards the Heat Equation 

In this section we will show that the flow of the aggregation-diffusion equation (1.1) behaves like the heat 
equation for large times provided we are under conditions for which the L^ decay of Theorem 3.1 holds. 
We prove it by two different arguments; though the results obtained by following these two strategies 
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are similar, we have kept both of them because they give a better understanding of the behavior of the 
equation. The proof in section 4.1 is based on direct estimates of the bihnear form (2.12). It is quite 
straightforward and gives a better result in the one-dimensional (N = 1) case. On the other hand, the 
argument given in section 4.2 is based on the self-similar change of variables and entropy arguments as 
in [17, 14], which we briefly describe now. 

We consider the standard self-similar scaling [16, 17] as is usually done for the heat equation to pass 
to the corresponding Fokker-Planck equation: 



/(s,y) = e%(^(e2^-l),e^y), (4.1) 



or, equivalently. 



p{t, x) = {l + 2t)-^/^f Q log(l + 2t), (1 + 2i)-^/'x) . (4.2) 
Then / satisfies the equivalent rescaled equation 

dsf = V, • (yf) + Ayf + Vy ■ {fiVyW * /)), (4.3) 
with initial data /(0,y) = po{y) if and only if p is a solution to (1.1). Here, we write 

W{s,y) :=W{ye'), so VyW{s,y) := e'VWiye'). (4.4) 
We consider the entropy functional 

H[f]:= I (^/log/+M!/^ dy. (4.5) 

We show below that H[f] converges exponentially to 0. The convergence for p obtained in this way 
is stronger, and requires stronger conditions on the initial data. Also, observe that in self-similar vari- 
ables the nonlinear term transforms to the time-dependent term involving W, where the asymptotic 
simplification becomes apparent. 

The results in this section are gathered in the following theorem. We recall that G is the fundamental 
solution of the heat equation, defined in (2.5). 

Theorem 4.1. Let the interaction potential W & Ci L"^ such that VW G {L^ n L°°)^ . Consider p the 
solution to (l.la)-(l.lb) with data po £ H L°° n H^, po > 0, constructed in Theorem 2.8. We denote 
M := \\po\\- Assume also one of the smallness conditions in Theorem 3.1. 

i) Then there exists C > such that: 

(a) For N = 1, 

\\pt- MG{t)\\^<^^^ for all t> I. (4.6) 
vt 

(h) For N> 2, 

-MG(t) II ^ < forallt>l. (4.7) 
\/t 

a) Assume further that |x|^po £ L^(M^). Then there exists C > such that: 
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(a) For N = 1, 

H[f] < Ce-' for all s >0. (4.8) 

(b) For N = 2, 

H[f] <C{l + s) e-"^' for all s > 0. (4.9) 

(c) ForN> 3, 

H[f]<Ce-^' forallsX). (4.10) 

Remark 4.2. We notice that under the above hypotheses Theorems 2.9 and 3.1 apply to the solution p 
in the statement. 

Remark 4.3. It is well-known that (4.10) directly implies, through the Csiszdr-Kullback inequality [33, 34] 
and the change of variables (4.1), that 

\\pt- MG{t)\\^<Ct~^/'^ forallt>0. 

This gives the same order of convergence as (4.7) for any N >3. However, the corresponding results for 
N = I or 2 are weaker than (4.6). In summary, point ii) above, which uses entropy /entropy dissipation 
tools, gives a worse rate of convergence than point i), whose proof uses a direct argument. However, one 
has to take into account that point ii) bounds a nonlinear quantity, the logarithmic entropy, essentially 
a quadratic functional at first order in expansion around the Maxwellian. Therefore, it is not surprising 
that we get a weaker result by estimating a quadratic functional. In fact, a direct argument in the spirit 
of point i) estimating the difference in Lp' of the solution with respect to the fundamental solution of the 
heat equation leads to similar rates of convergence as the entropy argument. 

4.1 Direct argument 

Let us prove point i) oi Theorem 4.1. We first give a technical lemma that will be used below: 
Lemma 4.4. For a > 1, we have the following bound 

f\t-s)-^l\l + s)-ds<\ # for all t> I, 

Jo y log(t) ifa = l 

where C > is some number which depends only on a. 
Proof. 

ft i-t/2 ft 

/ {t - s)-^/\l + s)-'' ds = / (l + s)-"(is+ / {t-s)-^/\l + s)-''ds = I + II 

Jo Jo Jt/2 

To estimate I we proceed as follows, bounding (t — s)~^/^ by (t/2)~^/^: 

ft/2 ( _C_ 

1= / (t - sr'/\i + s)-'^ ds < { 

Jo [ ^T72log(t+l) 

The estimate of II follows by 

// = [\t- s)"^/^(l + s)"" ds<{l + t/2)- 

Jt/2 

Adding the estimates for I and II gives the conclusion of the lemma. □ 



if a > 1 
if a = 1 



for alH > 1, 



/(.-.)- 



^'^ds< 



C 

t^ 
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Proof of point i) of Theorem 4-1- Taking into account (2.6), we get ||VG(t)||i < CNt~^^'^, for some con- 
stant Cn > depending only on the dimension. Using this estimate on the second term in the right-hand 
side of the Duhamel formula (2.8), we get 



hit) :-- 







{vG{t - s, •) * (pis, ■)iVW * p){s, •))) (x) 



ds dx 



< / ||VG(t-s,-)||i \\pis,-)iVW*p){s,-)\\ids 
Jo 

< Cn fit - s)-^l^ \\p{s, OlblKVt^ * p){s, ■)h ds 

Jo 

<Cn [ (t-s)"'/' \\p{s,-)h\\W\\i\\Vp{s,-)hds 
Jo 

< c\\W\\i f\t - s)-'/' (1 + s)-W2+V2) ds, 

Jo 

where we have used Theorems 3.1 and 3.5 for m = 1. Now, using Lemma 4.4 for > 2, we obtain 

h{t) < C\\W\\it-^/^ {t > 1), 

for some constant C > which depends on the dimension A^. From (2.8), this gives 

\\pt - e*^po||i < C\\W\\it-^/^ for ah t > 1. 

Using the known asymptotic behavior of the heat equation, for instance [18], the claim of the Theorem 
follows. In the case A'^ = 1, the same reasoning using the second bound in Lemma 4.4 gives 



\pt - e*^po||i < C\\W\\it-^/^ logt for all t > 1, 



which proves the result. 



□ 



4.2 Entropy argument 

We prove point ii) of Theorem 4.1: 

Proof. In addition to the entropy (4.5), we define as usual the entropy dissipation of the linear Fokker- 
Planck equation as 



D[f] ■■= J f 



VI ^ + log/ 



dy. 



As it is classically known [15], the evolution of the free energy for the equation (4.3) can be obtained 



as 



d 
ds 



H[f] + - II W{s, X - y)f{x)f{y) dx dy 



f 



V [^+logf + W*f 



dy + - 
^ 2 



dsW{s, X - y)f{x)f{y) dx dy, 
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and expanding part of the square in the first part yields the term D[f] on the right-hand side of the 
equation 



ds 



H[f] + \n W{s, X - y)f{x)f{y) dx dy 



-D[f]- J f VW*f dy 

+ 2 y / V + log • (VW *f)dy + \jj dsW{s, X - y)f{x)f{y) dx dy 

=:-D[f]+T2 + Ti + T^. (4.11) 
Recall that due to the classical Logarithmic Sobolev inequality [32, 31], we have 

2H[f] < D[f]. (4.12) 
We show that terms other than —D[f] decay in time at least like e~*. For the term we have 

73 = 2^ fy{VW*f)dy + 2 J V f ■ {VW * f) dy =: T^i + T32, (4.13) 

We prove that both T31 and T32 decay exponentially: for T31, using the Cauchy-Schwarz and the Young 
inequalities. 



Tsi < 2 ||y/||2 \\W * V/II2 < 2 \\yf\\, \\Vf\\, \\Wt = 2Ce-^' \\yf\\, \\Vf\\, < Ce 



~Ns 



(4.14) 



Here, we used that \\W\\i = e~^^ ||W^|li due to (4.4). The last step follows since both ||y/||2 and ||V/||2 
are uniformly bounded for all times by Theorem 3.5, (3.14) and the change of variables (4.1)-(4.2). By a 
similar argument for T32, we have 



T32 = 2 / V/ • (ly * V/) < 2 WVfWl \\W\\i = 2Ce-^^ ||V/||2 < Ce- 



(4.15) 



Using that ||V/||2 is uniformly bounded in time due to Theorem 3.5 through the change of variables 
(4.1)-(4.2). Note that, T2 < 0. For r4 we use that dsW{s,y) = y ■ VyW{s,y), the Cauchy-Schwarz and 
Young inequalities as above, to get 



7^4 = f{dsW*f)dy = ^ I fy{VW*f)dy = ^Tsi<Ce-^' for s > 0. 
On the other hand, the integral on the left-hand side of (4.11) can be bounded as follows 



W{s, x - y)f{x)f{y) dx dy 



<\j f{\W\*f)dy<]^ 



l\\W\\i<Ce-''\ 



(4.16) 



(4.17) 



Hence from (4.11), using (4.12), (4.13), (4.14), (4.15), (4.16), and (4.17), we obtain for all s > that 



^ (H[f] + f{W*f) dy^ < -2 (^H[f] + ^ I f(W*f)dy]+ Ce"^^^ + j fiW * f) dy 

< -2 (H[f] + 1 I f(W*f)dy]+ Ce-^^ 
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We remark that we have kept the term | Jf W{s, x — y)f{x)f{y) dx dy inside the time derivative to avoid 
the appearance of the time derivative of /. Hence, from this differential inequahty we deduce that there 
is some constant C which depends only on the initial condition /o such that 

H[f] + \j f{W * f) < Ce"'^^^^^'^}^ for N^2, 

with C(l + s)e~'^'^ on the right-hand side for N = 2. This implies, using again (4.17), 

H[f] < Ce-'^''i^^'2}s ^ _^ 2, 

for some other number C depending only on /q. Again, the right hand side should be substituted by 
C(l + s)e~^^ for N = 2. The final steps are just classical implications of this entropy control. Denoting 
the Maxwellian by 

A^(y) = (2vr)"^/2exp (-M!^ 
By means of the Csiszar-Kullback inequality, see [31, 33, 34], we get 

\\f{s) - MMis)\\l < CiHifis)) - H{M{s))) < Ce— ''^W^ls^ 
from which the announced result follows through the change of variables (4.1)-(4.2). □ 
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